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Establishment of cell polarity is crucial for many biological processes including cell migration and 
asymmetric cell division. The establishment of cell polarity consists of two sequential processes: an 
external gradient is first sensed and then the resulting signal is amplified and maintained by 
intracellular signaling networks usually using positive feedback regulation. Generally, these two 
processes are intertwined and it is challenging to determine which proteins contribute to the 
sensing or amplification process, particularly in multicellular organisms. Here, we integrated 
phenomenological modeling with quantitative single-cell measurements to separate the sensing 
and amplification components of Wnt ligands and receptors during establishment of polarity of the 
Caenorhabditis elegans P cells. By systematically exploring how P-cell polarity is altered in Wnt 
ligand and receptor mutants, we inferred that ligands predominantly affect the sensing process, 
whereas receptors are needed for both sensing and amplification. This integrated approach is 
generally applicable to other systems and will facilitate decoupling of the different layers of signal 
sensing and amplification. 
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Introduction 

Cell polarity is important for many cellular functions, such as 
migration, axis formation, and asymmetric cell division 
(Strutt, 2001; Nelson, 2003; Silhankova and Korswagen, 
2007} . Conceptually, the estabhshment of cell polarity in most 
cells can be thought of as a two-step process: sensing followed 
by amphfication (Weiner, 2002). The sensing machinery is 
used to detect an external cue. Examples of external cues for 
unicellular organisms include a pheromone gradient that 
induces cell polarity in budding yeast and cAMP gradients that 
directs the chemotactic response of Dictyostelium discoideum 
(Klein et al, 1988; Arkowitz, 1999). Remarkably, many cells 
are able to polarize in response to very shallow chemoat- 
tractant gradients as small as a 1 % change in concentration 
across the cell diameter (Zigmond, 1977; Baier and 
Bonhoeffer, 1992). This shallow gradient induces a steep 
intracellular gradient of signaling molecules and cytoskeleton 
components allowing the cell to polarize (Parent et aZ, 1998; 
Servant et aU 2000} . After the sensing machinery detects the 
external cue, an amplification mechanism sets in to convert 
this spatial information into a stable polarity axis. These 
amphfication mechanisms are often based on positive 
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feedback regulation (Meinhardt, 1999; Weiner, 2002). For 
example, during the establishment of cell polarity in budding 
yeast, activated Cdc42 orients the actin cytoskeleton and 
directs the delivery of more Cdc42 to membrane sites with high 
concentrations of the protein (Pruyne and Bretscher, 2000; 
Wedhch-Soldner et aU 2004). These vesicles are thought to 
contain Cdc42 therefore reinforcing the polarity axis via 
positive feedback (Wedhch-Soldner et aZ, 2003} . Amphfication 
through positive feedback has also been reported to be 
involved in planar cell polarity in multicellular organisms 
(Tree et aU 2002). In this case, positive feedback serves to 
polarize a field of cells by amplifying differences between 
protein levels on adjacent cell surfaces. Amphfication ensures 
robustness and is responsible for the dramatic sensitivity of 
gradient sensing and the generation of spontaneous cell 
polarization. 

Many mathematical models have been proposed to model 
the processes of cell polarization (Wedhch-Soldner et aU 2003; 
Jilkine and Edelstein-Keshet, 2011). Although these models 
have yielded valuable insights into the underlying molecular 
mechanisms required for sensing and amphfication, they 
cannot be used to separate and quantify the contributions of 
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the different molecular components to these processes, 
particularly in the in vivo context of a multicellular organism. 
Furthermore, construction of such mechanistic models 
remains challenging and is only possible for a few well- 
studied system due to the lack of detailed knowledge. On the 
other hand, phenomenological models that describe the main 
features of a system through a few important parameters have 
been used successful in understanding the general features of 
many systems (Alon, 2007; Mallavarapu et al, 2009). Here, 
through the integration of quantitative cell polarity measure- 
ments with a phenomenological model, we separated the 
contributions of Wnt hgands and receptors to sensing 
and amplification during establishment of polarity in the 
Caenorhabditis elegans' P cells. 

Wnt signaling plays an important role in regulating cell 
polarity in many organisms and Wnts have been found 
to be expressed in a gradient (Yoshikawa et aU 2003; 
Harterink et aU 2011). In C. elegans, loss of Wnt signaling 
leads to cell migration and cell division defects (Thorpe et aU 
1997; Whangbo etaU 2000; Silhankova and Korswagen, 2007). 
Many of these asymmetric divisions are controlled by the Wnt/ 
P-catenin asymmetry pathway. In this pathway, Wnts first 
regulate the polarity of the mother cell through asymmetric 
protein localization (Mizumoto and Sawa, 2007a}. After 
division, the asymmetric protein localization leads to different 
nuclear levels of POP-1 and SYS-1 in daughter cells (Phillips 
et aU 2007; Mizumoto and Sawa, 2007b}. This mechanism 
converts the polarity information set up in the mother cell into 
differential gene expression in the daughters (Mizumoto and 
Sawa, 2007b}. Although major progress has been made in 
uncovering new Wnt pathways and components, the exact 
functions of ligands and receptors remain unknown. In many 
cases, it was found that loss of Ugands and receptors leads to 
different phenotypes (Sternberg and Horvitz, 1988; Herman 
and Horvitz, 1994}. This is unexpected for components acting 
in the same pathway, suggesting that ligands and receptors 
may have different functions. But there are also cases where 



loss of receptors phenotypically mimick the effects of ligand 
loss. For example, both loss of ligands and receptors cause loss 
of polarity in the EMS cell (Thorpe et aZ, 1997}. Intriguing, the 
phenotype is less severe in the ligand/receptor double mutant 
than the hgand mutant (Bei et aU 2002}. This suggests that in 
the absence of hgands, receptors may be responding to signals 
that interfere with setting up of the correct polarity. 

Here, we found that Wnt signaling also regulates the 
division of the P cells, which are six pairs of symmetrically 
placed cells, P1-P12, found at the left and right sides of the 
worm after hatching. During the first larval stage, P cells 
migrate to the ventral side and form a single row of cells 
(Sulston, 1976}. After migration, they divide to produce 
anterior neuronal blast cells, Pn.a, and posterior epithehal 
daughters, Pn.p, 71 = 1, 2,... 12 (Figure lA}. P cells provide a 
convenient model system to study cell polarity for the 
following reasons. First, the estabhshment of cell polarity 
leads to a clear observable phenotype of different cell fates 
(neuronal versus epithehal} in the daughter cells. Second, this 
is an approximate one-dimensional system with anterior- 
posterior polarity. Third, the presence of many P cells in the 
same animal, allow us to study how the estabhshment of cell 
polarity varies along the anterior-posterior axis. 

Results 

lag-2 and lin-31 reliably mark the Pn.a and Pn.p 
cells, respectively 

To identify a marker of P-cell polarity, we quantified the mRNA 
expression of a panel of 26 genes using single-cell transcript 
counting (Raj et aU 2008} . This panel includes genes that were 
previously reported to be expressed in P cells and their 
descendants and also genes from the major signaling pathways 
(Wnt, Notch, FGF, EGF, and TGF}. A set of about 48 single- 
stranded 20-mer oligonucleotides were designed for 
visualization of each transcript. These fluorescently labeled 
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Figure 1 Identification of Pn.a and Pn.p expression marl<ers. (A) Lineage diagram of P cells. P cells divide to form Pn.a and Pn.p. Pn.a divides to give rise to five 
different neurons in L1 . Pn.p does not divide till the L3 stage. (B) Plot of lag-2 (red) and lin-31 (green) mRNA in Pn.p versus Pn.a. (C) (Top, left) Image of L1 larva after 
in s/fu hybridization with lag-2 (re6) and lin-31 (green) probes. Dapi is shown in blue. (Top, bottom) Maximum projection image of multiple z-stacks that shows lag-2 (re6) 
and lin-31 (green) expression in Pn.a and Pn.p from the P3-P10 lineages. (Right) Zoom in image of a single stack showing individual lag-2 an6 lin-31 transcripts in P8.a 
and P8.p. 
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oligonucleotides are complementary to the transcript and bind 
each individual transcript. This becomes visible as a diffrac- 
tion-limited spot using fluorescence microscopy. Using a 
custom-written software, we manually segmented the indivi- 
dual cells and computationally determined the transcript 
number in each cell. The ratio of the expression in the Pn.a 
versus Pn.p daughters was used to quantify the specificity of 
the putative markers (Supplementary Table SI). We identified 
Zag-2, a transmembrane protein of the Delta/serrate/lag-2 
family, as the most promising Pn.a marker having the largest 
expression ratio (37.3). lin-31, a gene previously reported as a 
Pn.p marker (Tan et aU 1998), was found to be the most 
promising Pn.p marker with lowest expression ratio (0.08). 
Using these markers, the anterior and posterior daughters of 
the P3-P10 cells can be scored with high confidence (Figures 
IC and 2C; Supplementary Figure 1). 

Although both daughter cells initially inherit lag-2 mRNA 
from the P cell, lag-2 is rapidly degraded in the Pn.p cells 
(Figure 2C}. We quantified lag-2 and lin-31 expression by 
counting individual transcripts in the Pn.p and Pn.a cells 
excluding cells immediately after division. We observed that in 
wild-type (WT) animals lag-2 and lin-31 expression is 
mutually exclusive (Figure IB). In WT animals, the lag-2 
count in Pn.a cells is always much larger than the lag-2 count 
in Pn.p cells (Figure IB). 



Role of the Wnt/|3-catenin asymmetry pathway 
in regulating P cells' divisions 

Many asymmetric divisions in C. elegans are regulated by the 
Wnt/p-catenin asymmetry pathway through polarization of 
the mother cell before division (Mizumoto and Sawa, 2007b}. 
In this pathway, asymmetry localization of many regulatory 
proteins (including WRM-1 and APR-l) in the mother cell is set 
up by Wnt ligands and receptors, leading to different nuclear 
levels of POP-1 and SYS-1 in the daughter cells after division. A 
high level of POP-1 and low level of SYS-1 in the anterior 
daughter leads to inhibition of Wnt signaling, whereas a low 



level of POP-1 and high level of SYS-1 in the posterior daughter 
leads to transcription of Wnt responsive genes. The amount of 
nuclear POP-1 in the posterior daughter is regulated by MOM-4 
and LIT-1 proteins, which function in exporting POP-1 from the 
nuclear into the cytoplasm. 

To test if the Wnt/p-catenin asymmetry pathway plays a 
similar role in P cells, we examined the nuclear levels of POP-1 
in the daughter cells after division in a POP-GFP translational 
fusion strain. We found that the anterior daughter has a higher 
level of nuclear POP-1 than the posterior daughter (Figure 2C} . 
This differential level of POP-1 is also observed in other 
daughter cells regulated by the Wnt/(3-catenin asymmetry 
pathway (Mizumoto and Sawa, 2007b}. Next, we disrupted the 
POP-1 branch in the asymmetry pathway by a temperature 
shift in mom-4; lit-1 mutant worms (Takeshita and Sawa, 
2005). Similar to WT animals, both Pn.a and Pn.p in mom-4; 
lit-1 mutant worms inherited lag-2 mRNA from the P cells. But 
unlike WT animals, Pn.p cells continue to express lag-2 some 
time after division (Figures 2B and 3B). This lag-2 expression 
is sustained even after Pn.p cells go on to divide in LI. 
(Figure 2D}. In WT animals, Pn.p cells do not divide until the 
L3 stage. This suggests that upon inhibition of the Wnt/ 
(3-catenin asymmetry pathway, the Pn.p cells have adopted 
Pn.a-like fates where they expressed lag-2 and divide in LI. 
This effect is similar to that observed in other cells whereby 
disruption of the Wnt/p-catenin asymmetry pathway caused 
the posterior daughter to take on anterior cell fate (Bertrand 
and Robert, 2009; Gleason and Eisenmann, 2010). These 
experiments confirm the role of Wnt/ (3-catenin asymmetry 
pathway in polarizing P cells and show that lag-2 expression in 
the daughter cells after division is a good reflection of the 
polarity of the mother P cell before division. 



Notch signaling is not involved in setting up P-cell 
polarity 

Since Pn.a and Pn.p express high levels of lag-2 and lin-12, 
respectively (Figure 2A}, we tested for the role of Notch 
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Figure 2 Inhibition of Wnt/p-catenin asymmetry patliway. (A) In WT animals, lag-2 transcript (red) and lin-12 transcript (green) are found in the Pn.a and Pn.p 
respectively. (B) In mom-4; lit-1 worms that have been transferred to 25°C for 9h prior fixation to inhibit the Wnt/p-catenin asymmetry pathway, /a^^-^ transcript is 
observed in both the Pn.a and Pn.p. (C) POP-1 localization is high in the Pn.a and low in Pn.p after P-cell division. (D) Both the Pn.a and Pn.p divide in L1 , leading to the 
increase in number of cells in the ventral nerve cord. 
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signaling in P-cell polarity. To determine if Notch signaling 
may be required for maintaining different fates in the Pn.a and 
Pn.p, we inhibited Notch signaling using a lag-1 temperature- 
sensitive strain (Qiao etaU 1995) . However, we did not observe 
any aberrant development in the Pn.a and Pn.p cells. We also 
examined animals with lin-12 loss of function mutation, 
lin-12(n941), and semi-dominant mutation, lin-12(nl37), and 
did not observe any change in cell fates and the expression of 
lag-2 in Pn.a and Pn.p cells (Greenwald et aU 1983). These 
results suggest that notch signaling is not involved in setting 
up P-cell polarity. 



Mutations in Wnt ligands induce polarity reversal, 
whereas mutations in Wnt receptors induce 
polarity loss 

To explore which Wnt ligands and receptors are involved in 
P-cell polarity, we examined lag-2 and lin-31 expression in Pn.a 
and Pn.p cells in ligand and receptor mutants. WT polarity was 
observed in most single- and double-ligand mutants, consistent 
with previous observations on the redundant role of the Wnt 
hgands (Zinovyeva et aZ, 2008} (Supplementary Figure S2}. 
However, in the double-ligand mutant [egl-20; cwn-1), triple- 



ligand mutant iegl-20; cwn-1; cwn-2), and quintuple-ligand 
mutant iegl-20; cwn-1; cwn-2; lin-44; mom-2), some Pn.p cells 
expressed lag-2 and their corresponding Pn.a cells expressed 
lin-31, suggesting that the P cell had been polarized in the 
opposite direction (Figure 3C; Supplementary Figure S3}. We 
found that in the triple-hgand mutant, 47.9% correctly 
polarized, 44.0% exhibited a polarity reversal, and 8.1% 
showed similar levels of lag-2 in both daughters indicative of a 
symmetric P-cell division and therefore loss of P-cell polarity. A 
very different phenotype was observed in the Wnt receptor 
mutants (Figure 3D}. In the mom- 5; lin-17 mutant, lag-2 was 
often observed in both daughter cells (53.5%} and a smaller 
fraction of cells correctly polarized (38.6%} or reverse 
polarized (7.9%}. The different phenotypes observed in the 
ligand and receptor mutants suggest that they may play 
different roles in establishing cell polarity. To determine if the 
different phenotypes observed is due to the different roles of 
receptors and ligands in sensing and amplification and to 
further quantify the phenotypes presented in Figure 3, we 
constructed a phenomenological model. 

Cell-cell signaling has been found to be important in planar 
cell polarity models (Tree et aU 2002} . Since pairs of P cells 
(P3/P4, P5/P6, P7/P8, and P9/P10} are in contact before 
division, cell-cell signaling may play a role in their divisions. 
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Figure 3 lag-2 expression in different mutants. Plot of lag-2 mRNA counts in Pn.p versus Pn.a in WT animals. (A) In animals with inhibition of Wnt/p-catenin 
asymmetry pathway, mom-4; lit-1 (B), multiple-ligand mutant, egl-20; cwn-1; cwn-2, (C) and double-receptor mutant, mom-5; Hn-17(D). Dotted lines divide the plot of 
/a^f-^ expression in Pn.p versus Pn.a into regions for correct polarity, loss of polarity and reverse polarity. Fractions of cells that fall into these three regions are indicated 
on the plots. 
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If cell-cell signaling occurs between pair of P cells, we would 
expect the polarity of division between these pairs of P cells to be 
correlated. In WT animals, all the divisions are correct hence 
WT animals cannot be used for this analysis. Hence, we 
calculate correlation in the triple-ligand mutant [egl-ZO; cwn-1; 
cwn-2) where many of the divisions are polarized in the 
reversed direction. If there is cell-cell signaling, we would 
expect pairs of P cells to both divide correctly, leading to both 
Pn.a daughters having high lag-2 levels, or to both divide 
with reverse polarity, leading to low levels of lag-2 in both 
Pn.a daughters. When we examined the lag-2 expression levels 
in the anterior daughters for the anterior P cell versus the 
posterior P cell in each pair (Supplementary Figure 4), no 
significant correlation is found between their expression levels 
(correlation value = -0.029, P-value = 0.70). This suggests 
that cell- cell signaling does not play a significant role in 
P-cell division. 

A phenomenological model for cell polarity 

The polarization state of the P cell is denoted by the dynamic 
variable 9(0 , which is an angle that ranges between 0° (reverse 
polarity) and 90° (WT polarity). We assume that the P cell has 
initially no polarity (9(f= 0) =45°). For subsequent times, 
until the P cell divides at time T, we describe the dynamic 
behavior of 9(0 as a one-dimensional random walk. In the 
simplest model that only involves the sensing of an external 
gradient, the probability of taking a right step (increasing 9; 
towards WT polarity) or left step (decreasing 9; towards 



reverse polarity) is ((1 +g)/2)rTAt and ((1 -g)/2)rTAf, respec- 
tively. Here, g reflects the external gradient that ranges from 
- 1 to 1, r-r is the total reaction rate that sets the timescale of 
the dynamics, is the time step of the numerical simulation. 
To understand the behavior of the model, we ran stochastic 
simulations to determine the behavior of individual cells. The 
chemical master equation is also solved to obtain the 
distribution for a population of cells at time T. In absence of 
an external gradient (g = 0), 9(0 has equal probability of 
increasing or decreasing at each step of the simulation 
(Figure 4A) and performs an unbiased random walk resulting 
in an approximately Gaussian distribution (symmetric about 
45°) at the time of P-cell division (Figure 4B, green). If g>0 or 
g<0, the probabilities of taking a left and right step are 
unequal, resulting in a biased random walk towards 90° (blue) 
or 0° (red), respectively (Figure 4A and B). To summarize the 
information encoded in the histograms of 9(7) for each set of g 
and r-r, we determine the fraction of the cells that undergoes 
correct, reverse or loss of polarity based on the following 
classification rules: 9(7) > 75° (correct polarity), 9(7) < 15° 
(reverse polarity), and 15°<9(7)<75° (loss of polarity). A 
coordinate on the phase diagram corresponding to the 
particular set of g and r^, will be colored with different 
intensities of blue, red, and green depending on the fractions of 
correct, reverse, and loss of polarity, respectively. (Figure 4C) 
As expected for g>0 and high r^, values of most cells are 
correctly polarized. For values of g close to 0 and low values of 
Tt, most of the cells are unpolarized. However, this simple 
model is unable to explain the coexistence of correct and 
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Figure 4 Stochastic modeling. Without feedbacl< (A) Schematics showing the dependence of and rp on g. For = 0, = Tp and for 0 < 1 , > Tp. Values of 
and Tp do not depend on 6a. (B) (Left) Stochastic simulation of 6a from f = 0 to f = T (division of P cell) (Right) Histograms of 6a(7) for constant rT= 0.05 with different 
values of 9^ = 0.3 (blue), g = Qi (green), and 5^= - 0.3 (red). (C) Phase diagram showing the parameters space of g^and r^giving rise to correct polarity (blue/C), reverse 
polarity (red/R), and unpolarized divisions (green/U). Dotted lines show the region of parameter space giving rise to 35% correct polarity (top) or 35% reverse polarity 
(bottom). With feedback (D) Schematics showing the dependence of and rp on g and 6. For g = 0, ra = rp at 63 = 45°. For 0 <g'< 1 , ra>rp at 6a = 45°. For both 
cases, Ta increases with 6a whereas rp decreases with 6a. (E) (Left) Stochastic simulation of 6a from t=0 to t=J (division of P cell) (Right) Histograms of 6a(T) for 
5^ = 0.3 and rf= 0.06 (blue), for 9^= - 0.3 and r^= 0.06 (red), \or g=0 and rf= 0.06 (magneta) and \or g = 0 and rf= 0.01 (green). (F) Phase diagram showing the 
parameters space of g^and giving rise to correct polarity (blue/C), reversed polarity (red/R), unpolarized divisions (green/U), and coexistence of correct and reverse 
polarity (magneta/C/R). Dotted lines show the region of parameter space giving rise to 35% correct divisions (top) or 35% polarity flips (bottom). The intersection of the 
two dotted lines shows the coexistence of correct and reverse polarity. 
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reverse polarity observed in the triple-ligand mutant 
(Figure 3C}. 

To expand beyond this simple model, we included a 
mechanism that reinforces a deviation from the non-polarized 
state. Biologically, this reinforcement could be estabhshed by, 
for example, positive feedback regulation. To introduce 
amplification, we let the probability of taking a right step 
increase with 9 and the probability of taking a left step decrease 
with 9. The simplest way to reinforce deviation is to introduce 
linear dependencies on 9 into the probabilities of taking a right 
or left step. In this expanded model, the probability of taking a 
right step (increasing 9; towards WT polarity) or left step 
(decreasing 9; towards reversed polarity) is ((1 +g)/2)(9r^+ 
rij)At and ((1 -g)/2)([90° - 9]r/-+ r5)Af, respectively, g per- 
forms the same role here as in the simple model and 
characterized the gradient. On the other hand, r-r is replaced 
by two rates r^^and r^. In the absence of a gradient (g = 0), 9(0 
has equal probability of increasing or decreasing when the cell 
is unpolarized (9 = 45°) similar to the behavior of the simple 
model. However, if there is any fluctuation that drives the cell 
away from the unpolarized state, this fluctuation will be 
reinforced. In other words, the probability that 9 will move in 
the same direction as the fluctuation does, is larger than the 
probability to move in the opposite direction of the fluctuation 
(Figure 4D and E, green traces) . The parameter Vf quantifies 
this reinforcement. A basal rate independent of 9(0, r^, is 
included and set to be a constant. The external gradient is 
superposed on this process and further biases the random 
walk (Figure 4D and E, blue, red, and magenta traces). 
Another way to view the polarization process is that initially 
when the cell is unpolarized (9 = 45°), the contributions of r^to 
the probabilities of taking a right or left step are equal. Hence, 
the sensing process, represented by g, plays a more important 
role is setting up polarization. Once sensing has occurred and 9 
is no longer close to 45°, contributions of amphfication 
represented by will become important and act synergistically 
with the gradient to further set up the polarization. A model 
assuming additive effects of gradient sensing and ampUfica- 
tion is also unable to produce the coexistence of correct and 
reverse polarity (Supplementary Figure S5) . 

We repeated the analysis for the expanded model, setting 
as a constant and found a region in parameter space (g versus 
Tf), which allows for the coexistence for correct and reverse 
polarity (Figure 4F, magenta), which was absent in the simple 
model (Figure 4C) . We then used this model to further quantify 
the experimental results and explore where the different 
mutants are located in this parameter space. Note that the 
present model is a purely phenomenological model (in 
contrast to a mechanistic molecular model) that allows us to 
extract parameters, g and r^, from the experimental data. 



Loss of Wnt ligands reduces the parameter g, 
whereas loss of Wnt receptors results in a 
decreased value of both g and 

To compare the predictions of the model to the experimental 
data, we converted the lag-2 transcript count in each pair of 
daughter P cell to a single angle Qiag-z (Figure 5A). From the 
earher experiments, we observed that inhibition of the Wnt/p- 

6 Molecular Systems Biology 2013 



catenin asymmetry pathway, which sets up polarity of the 
P cell, led to expression of lag-2 in both daughter cells. Hence, 
diag-2 is a good approximation of the polarization state of the 
mother P cell at the time of division 9(7) as calculated by the 
model above. 

To determine the value of r^, we examined the distribution of 
diag-2 in the mom-4; lit-1 strain and observed a Gaussian 
distribution (Figure 5B). Although ligands and receptors are 
functioning properly in the mom-4; lit-1 strain, Pn.a and Pn.p 
are unable to execute different cell fates as the MOM-4 and 
LIT-1 proteins, required for the daughter cells to have different 
amount of POP-1, lost their function. We observed that the 
distribution for diag-z is centered at 60°. The bias of the diag-z 
distribution is likely to be due to the differential amount of 
SYS-1 in the two daughter cells. Since information determining 
through sensing and amphfication are not conveyed effectively 
to the daughter cells as the POP-1 branch is inhibited, we can 
use this mutant to obtain an estimate for = 0.006 by setting 
rf=0. This is likely to be an overestimation as the SYS-1 
branch is unaffected. Using this estimation for the value of r^,, 
we determined the values of g and Vf that yields maximum 
likelihood fits for the WTand mutants strains' distributions. 

When we fit the WT distributions, we obtained large 
uncertainties in the fit parameters. This is expected because 
the parameter range over which WT polarity is observed is 
large. Although the WT distribution cannot be fitted uniquely, 
we know that the parameters lie in the parameter space where 
100 % correct divisions are observed which is bounded by the 
dashed line in Figure 5D and E. Next, the experimental 
distribution of diag-z for the different P cells in the ligand and 
receptor mutants (Figure 5C, blue histograms) were fit to the 
model (Figure 5C, red lines). The distributions for each of the 
individual P cells are significantly different; hence, the fit 
parameters g and r^^ were determined for each P cell in the ligand 
mutants and receptor mutants independently. The model is 
able to reproduce the main features of the experimental 
distributions. Interestingly, these two parameters segregated 
out in parameter space (Figure 5D). We found that the gradient 
parameter g was always lower compared with WT (Figure 5D), 
suggesting that both Wnt ligands and receptors are important 
for the sensing of the gradient. However, the feedback 
parameter Vf was significantly smaller when receptors were 
mutated compared with ligand mutations, suggesting that the 
receptors might be involved in amplifying the signal. 



Decreasing ligand results in symmetric divisions 
at low rf 

From the above analysis, we inferred that Wnt hgands are 
primarily involved in sensing whereas Wnt receptors function 
in both sensing and amphfication. Mutations in Wnt hgands 
primarily lead to polarity reversals whereas mutations in Wnt 
receptors cause loss of polarity. Very similar results have been 
reported for many tail blast cells in C. elegans, including the 
T cell, suggesting that our model maybe generally apphcable 
(Sternberg and Horvitz, 1988; Herman and Horvitz, 1994). 
However, in contrast, it has been reported for the EMS cell in 
C. elegans that loss of ligands causes loss of polarity apparently 
inconsistent with our interpretation (Thorpe et aZ, 1997). 
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However, upon closer inspection of the polarity phase diagram 
(Figure 5E), we hypothesized that in the EMS cell, the WT 
parameters are likely to be different. For example, the EMS cell 
might experience a different gradient or has reduced amplifi- 
cation potential compared with the P cells. 

The parameter range over which WT polarity is observed is 
large and is bounded by the dashed line in Figure 5D and E. If 
we would reduce Vf slightly so that most P cells still correctly 
polarize (moving from WT to C in Figure 5E} and subsequently 
remove ligand (moving from C to U in Figure 5E} it should be 
possible to induce loss of polarity upon loss of Wnt ligand. To 
test this prediction and access the parameter space of low r^, 
we looked for mutants for which Vf is lower than WT but 
sufficiently high to yield correct divisions. We chose the mom- 
5 single receptor mutant as our candidate as there are many P 
cells with correct divisions in this mutant. In the mom-5 
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mutant, divisions of the P3 and P4 cells are symmetric whereas 
the divisions of the P5-P10 cells are correct (Figure 6A}. 
Reducing hgand levels in the mom-5; cwn-1 strain led to 
symmetric divisions in the P5 and P6 cells (Figure 6A} 
confirming our prediction. This behavior was also reproduced 
in the mom-5; egl-20 and mom-5; egl-20; cwn-1 mutants, where 
divisions in all P3-P10 cells were affected (Figure 6A}. 



Positive correlation between receptor levels and tf 

In the mom-5; egl-20 and mom-5; egl-20; cwn-1 strains, where 
divisions in all P3-P10 cells were affected, we observed a high 
fraction of symmetric divisions in the P3-P6 cells, lower 
fractions of symmetric divisions in the P7-P8 cells and 
coexistence of correct and reverse polarity in the P9 and PIO 
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cells (Figure 6A}. It is intriguing that all these different 
phenotypes were observed in the same mutant strain and even 
in the same mutant animal, suggesting that the values of Vf 
may vary significantly among P cells in the same animal. We 
hypothesized that the differences in Vf are due to the different 
expression levels of the receptors. We measured transcrip- 
tional levels of the receptors genes, mom-5 and lin-1 7, using 
single molecule FISH in the P cells and found that anterior P 
cells expressed higher levels of mom-5 whereas posterior cells 
expressed higher levels of lin-17 (Figure 6B}. We found a 
positive correlation between the lin-17 count and the 
parameter Vf (Figure 6C) for both the mom-5; egl-20 and 
mom-5; egl-20; cwn-1 strains. Similarly, a positive correlation 
was observed between the mom-5 mRNA level and Vf in the 
lin-17; egl-20; cwn-1 strain (Figure 6D). These positive 
correlations obtained between receptor levels and Vf provide 
quantitative support for our earher conclusion that receptors 
may play an important role in amphfication. Furthermore, 
since differences in receptors level could explain most of the 
differences in r^, it also demonstrates that the ligands' 
contribution to amphfication is less significant. 

Discussion 

By combining quantitative single molecule transcript counting 
with phenomenological modeling, we studied the effects 
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of hgand and receptor loss on P cells' division. We found that 
hgand affect primarily the sensing process whereby loss of 
hgands lead to polarity reversals. Single molecule FISH 
experiments on Wnt ligands show that ligands are expressed 
in a gradient, cwn-1 and egl-20 are higher in the posterior 
region whereas cwn-2 is higher in the anterior region 
(Harterink et aU 2011). These ligand gradients could poten- 
tially serve as the gradients for polarizing the P cells. Indeed, 
Wnt hgands has been shown to act instructively (Goldstein 
et aU 2006). But there are also cases where hgands act 
permissively (Whangbo et aU 2000) . To determine if hgands 
act instructively or permissively, we expressed cwn-2 poster- 
iorly using an egl-20 promoter or uniformly using a myo-3 
promoter in a cwn-1; cwn-2; egl-20 triple-ligand mutant 
(Supplementary Figure 6) . We found that ubiquitous expres- 
sion of cwn-2 using both promoters are able to rescue the 
polarity reversals, suggesting that it is the ligand level rather 
than the gradient that is important for the sensing process. 
However, the recent identification of Wnt inhibitor, sfrp-1, in 
the anterior region of the worm suggests that hgand profile is 
not shaped by ligand expression alone (Harterink et aZ, 2011). 
Anterior expression of sfrp-1 could potentially convert uniform 
hgand expression into a hgand gradient. Identification of more 
Wnt components and quantitative measurement of ligand 
protein profiles wih be necessary to resolve this issue. 

Unlike the hgand mutants, polarity loss is observed in the 
receptor mutants. We reconcile the different phenotypes 
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observed in ligand versus receptor mutants. Ligands 
affect primarily the sensing process whereas receptors 
are needed for both sensing and amphfication. This suggests 
that receptors can signal to the downstream proteins needed 
for amphfication regardless of whether they are Ugand-bound 
whereas the function of hgand is to convey information 
about the directionality of polarization to the cell through 
the receptors. One possible molecular mechanism is that 
initially equal amount of receptors bind to both the anterior 
and posterior sides of a cells. The presence of a hgand 
gradient would bias the amount of receptors on one side 
of the cell through mechanisms like reduction in degradation 
rates for ligand-bound receptors. The higher amount of 
receptors on one side of the cell may lead to positive feedback 
to bring even more receptors to the same site. This would 
eventually cause the cell to become polarized and divide 
asymmetrically. In the absence of ligand, positive feedback 
could still lead to a higher amount of receptors on one side of 
the cell due to initial stochastic differences in the amount of 
receptors at each side of the cell. However in this case, some 
fraction of the cells will be polarized in the opposite direction, 
leading to existence of divisions of both correct and reverse 
polarity. When receptors are absent, there will not be positive 
feedback and the cells would divide symmetrically. The 
observation that receptors cluster asymmetrically before 
division is consistent with this proposed model (Wu and 
Herman, 2007). Future work is necessary to elucidate the 
molecular mechanisms for the differences in roles of Ugands 
and receptors. 

In our mom-5;wnt strains (Figure 6A}, we observed loss of 
polarity and polarity reversals for P cells with low and 
high Tf, respectively. This is reminiscent of the different 
phenotypes observed upon ligand loss in C. elegans. Loss of 
ligand leads to polarity reversals for many tail blast cells, 
including the T cell, and loss of polarity in the EMS cell. Our 
result suggests that these differences are due to low amplifica- 
tion strength in the EMS and high amplification strength in the 
T cell. It is interesting why some cells have high amplification 
strength whereas others have low amplification strength. We 
suggest that low amplification is sufficient for correct division in 
the EMS as the single MOM-2/ligand source is adjacent to the 
EMS and likely to provide a steep gradient across the cell. 
However, for P cells during post-embryonic development, 
the gradient is shaped by many different ligands and by 
many ligand- expressing cells all along the worm (Harterink 
et al 2011). Therefore, the gradient experienced by the 
P cells is likely to be shallower and more complicated to 
interpret. In this case, it will be beneficial for the cell to have 
higher amplification is to ensure robustness. Hence, we 
speculate that depending on the gradient and the complexity 
of the ligand sources, cells may express different receptor levels 
to balance the tradeoff between robustness of polarity and 
energy involved in high expression. 

Here, we have proposed a general framework for determin- 
ing the contributions of proteins to sensing and amphfication. 
This method overcomes the need for detailed knowledge 
of the biochemical reactions and network topology involved 
in many polarizing system and allows us to extract 
parameters that quantify the key processes of sensing and 
amphfication. 
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Materials and methods 

C. elegans strains and culturing 

General methods for culture, manipulation, and genetics of C. elegans 
were as described previously (Lewis and Fleming, 1995). Mutations 
used in this study were: LGI, lin-44[nl792) (Herman and Horvitz, 
1994), mom-5[or57) (Thorpe et al, 1997) in mom-5; lin-17 receptor 
strain and mom-5[gk812] for all other mom-5 strains; LGII, cwn- 
l[ok546) (Zinovyeva and Forrester, 2005), mom-4[nel539) (Takeshita 
and Sawa, 2005); LGIII, lit-l[tl512] (Takeshita and Sawa, 2005); LGIV, 
cwn-2[ok895] (Zinovyeva and Forrester, 2005), egl-20in585] 
(Zinovyeva et al, 2008); LGV, mom-2(ne874ts). Strains were cultured 
at 20°C except for the mom-4; lit-1 strain, which was cultured at 15°C 
and transferred to 25°C for 7 or 9h before fixation and the Wnt 
quintuple mutant, which was grown at 15°C and shifted to 25°C for 7 h 
before fixation. 



Single molecule mRNA FISH 

Probe design and hybridization to perform FISH for single transcript 
measurement in C. elegans larvae was performed as previously 
described (Raj et al, 2008). Animals were collected by washing plates 
with M9 and were fixed in 4% formaldehyde in 1 x PBS for 45 min. 
Fixed animals were permeabilized in 70% ethanol overnight. All 
probes for hybridization were coupled to either Cy5 (GE Amersham) 
or Alexa594 (Invitrogen), depending on the desired gene com- 
binations for image acquisition. Images were taken with a Nikon Ti-E 
inverted fluorescence microscope equipped with a 100 x oil- 
immersion objective and a Photometries Pixis 1024 CCD camera 
using MetaMorph software (Molecular Devices, Downington, PA, 
USA). Three-dimensional positions of bright fluorescent spots 
in each animal were detected with the aid of a custom program 
written in MATLAB, as described (Raj et al, 2008), which was later 
manually corrected for further accuracy. Nuclei were visualized 
with DAPI. 



Simulation of histograms of B 

The chemical master equation was solved to obtain the exact solution 
for each set of parameters. In the simulation, 9 takes values between 0 
and 90° and increases or decreases in steps of 1°. Next, a transition 
matrix, M, describing the rates of going from between the different 
values of 9 is constructed. Since there are 91 different values that 9 can 
take, the transition matrix is of size 91 x 91 . The rate of going from 9i to 
92 is represented by M(92 + 1, 9i +1) in the matrix. Furthermore, since 
9 can only increase or decrease in steps of 1°, only the diagonal terms 
and values next to the diagonal terms are non-zero in the transition 
matrix. 



For the simple model 

M(0 + 2, 0 + 1)= ^^Tt 

for 9 = 0°-89° 

M(0, 0 + l)=i^rT 

for 9 = r-90° 

M{e, 6)^ -(M(0-i, e)+M{e + i, e)) 

for 9 = 2°-90° 

The boundary conditions at 9 = 0° and 90° are reflecting. 

M(l, 1)= -M(2, 1) 

M(91, 91)= -M(90, 91) 
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For the expanded model 

M(0 + 2, e+l)=^^{9rf + rb) 

fore = 0°-89° 

M(6', e + l) = ^^i[90-e]rf + rt) 

fore = r-9o° 

M{e, e)= -{M{e-i, e)+M{e + i, e)) 

fore = 2°-90° 

Since 9 cannot takes values less than 0° and greater than 90°, the 
boundary conditions at 9 = 0° and 90° are reflecting. 

M(l, 1)= -M(2, 1) 

M(91, 91)= -M(90, 91) 

The distribution at time t is presented by and is generated by the 
following equation: 

Pt+i = il-MAt)Pt 

where I is the identity matrix. We start the simulation with Pf=o 
(9 = 45°) = 1 as the cells are initially unpolarized. The total simulation 
time T is underdetermined and is set constant to 1 for all the 
simulations and different mutants analyzed. At is the time step of the 
numerical simulation and has been chosen to be 10"'^ in our 
simulation. The resulting histograms are independent of the choice 
of At as long as At is small enough. 

Classification rules for correct polarity, reverse 
polarity, and loss of polarity 

The classification rules are 9(7) > 75° (correct polarity), 9(7) < 15° 
(reverse polarity), and 15° ^9(7) ^75° (loss of polarity). These rules 
are obtained by finding the angle that best separate the lag-2 
expression of Pn.a versus Pn.p between WT animals and the mom-4; 
lit-1 strain at 25°C (Figure 3 A and B). 



Stochastic simulation of single B(f) traces 

Monte Carlo simulation were performed to generate many single 9(f) 
trajectories. The initial condition 9(f=0) =45° was used as the cells 
are initially unpolarized. For every time t, a random number n is 
drawn. 9(f + At) will be updated based on the following conditions: 



Simple model 

e{t+At) = e{t) + i 

for 

1 +^ 
n< -^rjAt 

e{t+At)^e{t)-i 

for 

1 +g 1 +g 1 — g 

rjAt <n< rjAt H ^ rjAt 

e{t+At)^d{t) 

for 

1 +g . 1 -g . 

n> -r-iAt^ ^TtA^ 

2 2 



Expanded model 

e{t+At) = e{t) + i 

for 

n< ^-^{erf + r^)At 

e{t+At)^e{t)-i 

for 

^{evf + r,)At <n<^^{9rf + r,)At 
+ ^ {[90 -e]rf + r,) At 

d{t + At)^d{t) 

for 

n> ii^(6)r^ + rb)At+i^([90-6']r^ + r5)At 
The boundary conditions at 9 = 0° and 90° are reflecting. 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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